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A. Introduction 

The terminal part of the e"*" track (the positron blob) is formed during ionization slowing 
down and subsequent ion-electron recombinations produced by a positron. It releases up to 
1 keV of energy, which is converted into heat within few picoseconds. If a bulk temperature 
of a medium is below, but close enough to its melting point, some region of a substance may 
melt, yielding a peculiar temperature dependence of the lifetime (LT) spectra. 

We have estimated properties of the molten region with a help of macroscopic heat 
conduction equation and suggested a model describing temperature dependence of the ortho- 
positronium lifetime in frozen methanol, ethanol, butanol and water close to their melting 
points. 

B. General equations. Homogeneous medium 

In a homogeneous medium the heat conductivity equation for the e^ blob looks as follows: 

CpP^^^ = div(AVT) + g+(r, t), T{r, t = 0) = T^uik- (1) 

Here T{r,t) is the local temperature, Tbuik is the bulk temperature of the medium, Cp is 
its specific heat capacity, p is the density, A is the thermal conductivity. The second term 
in the RHS describes energy release by the positron when it creates the blob: g+(r, t) ^ 
WbiobG{r,a)f(t,T), where Wuob ~ 1 keV is the blob formation energy, G{r,a) = ^ 3/2 3 
describes spatial distribution of the released energy (a is about tens of A) and /(t, r) its 
temporal distribution (here r is the typical time of ion-electron recombination and should 
be of the order of 1-10 ps). We adopted that /(t, r) = exp ( — ^ ~2r^ ) /(v^^ttt), where 
r = 0.3 ps. Since a system is simultaneously solid and liquid we used the following method 



fl. 



for obtaining a numerical solution of Eq. ([T]) [ij. An additional contribution to c^iT^ was 
used to simulate a presence of the latent heat of melting (gm)- It is non-zero only in a small 
temperature interval AT around T^ (the melting point temperature) and its integral over 
temperature in AT range should be equal to g^- The value of AT (phase transition width) 
was arbitrarily fixed to 0.25 K. An additional contribution to Cp related with the latent heat 
was simulated by a Gaussian T-dependent function. Dependencies of c^iT^ as well as A(T) 
and piT^ are described by smooth functions over T as proposed in pL,]. Thermodynamical 



properties of methanol, ethanol, butanol and water used in present work are taken from 



C. Estimated properties of molten region 



Temperature profiles T{r, t) were calculated numerically for < r < 200 A and < t < 1 
ns. Eq. ([1]) was solved as 1-D problem with boundary conditions T{r = 200A) = T^uik with a 
help of PDEPE solver from Matlab. The phase of medium is deduced from the temperature of 
each spatial point, which can be higher than T^ + AT (liquid phase) or lower than T^ — AT 
(solid phase). 

Fig. [1] displays the maximum radius of the molten region, Rmax, vs. Tbuik and the lifetime 
tmax of the molten region (at t > tmax temperature of an any point of the medium is below 
T^ + AT). 
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FIG. 1: Lifetime of the molten region tm (left) and its radius Rm (right) for methanol, etanol, 
buthanol and water are displayed as a function of Tj, — Tm- 



D. Experimental data. 



Positron annihilation lifetime spectra in alcohols and in water were measured in Maria 
Curie-Sklodowska University. Each lifetime spectrum was deconvoluted into 3 exponents by 
means of LT 92 program. The data were processed at a fixed ratio of ortho-Ps to para-Ps 
intensities equal to 3:1. Fig. [2]displays temperature variation of the lifetime of the long-lived 
component of LT spectra vs. T close to the melting points. 

For alcohols one can observe a noticeable temperature shift between the beginning of 
r3(T) plateau and the melting point. For description of r^iT) dependencies the following 




FIG. 2: Experimental data for the lifetime of ortho-positronium component in alcohols and water. 

Temperature T of investigated medium is plotted relatively to Tm of certain substance. 
model is proposed. 



E. Theoretical model. 



Let us divide a space-time evolution of the Ps into two subsequent stages, namely, its 
formation and annihilation. The first stage is initiated by a heat generation during e^ 
ionization slowing down, when it forms a terminal blob. As a result, intrablob region may 
melt, so therein formation of the Ps bubble is possible. Typical duration of the first stage 
may be estimated as about tmax ~ 10 ps, so annihilation may be neglected. 

When temperature of the medium approaches to the melting point, concentration of free 
volume elements (FVE; or vacancy-like defects) rise extensively. These FVE trap Ps and 
yield its typical lifetime as about r^ ^ 1.2-1.5 ps. 

Frozen phase usually consists of domains which are disordered in a different extent (local 
density fluctuations, intercrystalline boundaries and so on). Local inhomogeneity of the 
medium leads to smearing of its melting temperature. It is the reason of the premelting 
effect (an appearance of some molten domains when Th^ik is slightly below T^)- Hence the 
formation of the Ps bubble state may occur in these premelted regions. 

Let us consider Ps formation during the first stage in terms of concentrations of different 
Ps states, namely, mobile quasi-free positronium qf-Ps, Cg/(r, t), Ps trapped in FVE of a 
sohd phase, c„, and Ps bubble state, c^, formed in the molten region of the blob and in the 
premelted domains. Evolution of these states is described by the following set of kinetic 



equations: 

0^ = D/^C^f - [\v + \b)Cqf, ^^ = \Cqf, — = XbCqf. (2) 

Here D ~ 0.03 cm^/sec is the diffusion coefficient of qf-Ps (we assume that the mean free 
path of Ps transport is about its de Broghe wave length and qf-Ps velocity is the thermal 
one), XviT) and A(,(T) are Ps trapping rates by FVE and premelted regions. 

At t = we have qf-Ps only. We adopt that its initial distribution is approximately 
Gaussuan and that it remains Gaussian later at any t, namely, Cqf{r,t) = nqf{t)G{r,t,D), 
where G{r,t,D) = '''^[f(^l/4g^^3/? • Hence nqf{t) = Pqf_Pse-(^''+^'')*, where Pqf_Ps is the 
formation probability of qf-Ps. Distributions c^ and c^ can be found from ([2]) and are equal 
to c^,b{r,t) = /q dt' X^^bCqf{r,t'). 

The first stage is terminated by a formation of the molten region with the radius, attaining 
its maximal value Rm{Tbuik) at tmax- By this time the fraction /^ of Ps atoms in premelted 
regions (i.e. in a hquid phase) is fm = X-^p {cqf{r,tmax) + c^{r,tmax) + Cb{r,tmax))d^ 



r. 



By the time tmax ~ 10 ps all these Ps states turn out to be in a bubble state (Ps bubble 



Is 



formation time is about 10 ps as well p]). After subsequent freezing of the molten region 
we assume that the bubble state does not collapse. Cavity survives after freezing with the 
Ps inside, so ortho-Ps lifetime remains the same. 

We assume that all qf-Ps atoms outside the molten region, remaining quasi-free by the 
time tmax, become trapped shortly by FVE or premelted regions (n„ or Ub) in proportion to 
Ai>/(At, + Afe) and Xb/{Xv + Xb). Therefore, total formation probability Pps of the Ps bubble 
state arises not only from the molten central region of the blob, but also from premelted 



domains, so that Pps = fm + J^^r Cb{r, tmax) + ATAr^'?/(^' ^ 



A„+A5^9/V'''™«^^ 



rf3, 



Annihilation of the free positrons and Ps states occurs on the second stage at t > tmax ~ 0. 
We assume that formation of para- and ortho-Ps spin states takes place in a conventional 
1:3 proportion (possible influence of intratrack chemical reactions which may change this 
ratio is neglected). Let ripPs (oPs) be the fraction of para-Ps (ortho-Ps) formed in a bubble 
state. Following equations allow to calculate a shape of LT spectrum according to proposed 
model: 

n+ = -A+n+, ra+(0) = 1 - Pqf.p^; rioPs = -Apo^oPs, ^oPs(O) = 3Pps/4; 

npPs = -ApPsTlpPs, '^pPs(O) = ^ + - / \C^{r,tmax) + \ ^^Cqf{r,tmax)] d^T] 
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FIG. 3: Temperature dependencies of trapping rates A„ (solid lines) and A^ (dashed lines) for 
methanol, ethanol, butanol and water. Temperature is given with respect to melting point of each 
substance. 
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These equations yield a four-component LT spectrum, where A4. ~ 2 ns^^ is the free 
e"*" annihilation rate, ApPs is the annihilation rate of para-Ps (we do not distinguish p-Ps 
either in a bubble state or not), Xpo is the pick-off annihilation rate of Ps in a bubble state 
and Xs is the o-Ps annihilation rate in free volume elements of frozen regions. However, by 
now experimental data were deconvoluted into three exponents. So we used the following 
expression to describe the long-lived component T3 of the spectra: 



Ts ^,;(0) + noPs(O) ^°n^(0) + noPs(O) ' 



(3) 



Implicitly two temperature dependent quantities, A^, and A^, enter this expression. It is 
reasonable to describe them by the Arrhenius-like law A^,,f,(T) ~ exp{—Ey^i)/T) and extract 
corresponding activation energies for free volume elements E^ and premelted regions Ef,. 

Fitting the model (|3]) to the experimental data (Fig. |2]) in the range close to melting 
points, we obtain temperature dependencies of trapping rates. Fig. [31 To eliminate corre- 
lation of parameters in the Arrhenius law the following procedure was applied. Melting is 
assumed to occur when the number of premelted regions (let a typical size of these regions is 
Rs) allows to fill the whole space, i.e. Xb(T = Tm)/ {inDRs] ■ ^t^R^s — ^- From this relation 
it is possible to get an additional constraint for Ab(T = T^). Furthermore, if we consider 
the energy Ei, as an energy necessary for melting of the media inside a premelted domain 



(|7ri?|), then Ef, = g^' |^-R|' Ps- Thus obtained relations for E;, and A{,(T = T^) are depen- 
dent on the only parameter Rs which has a meaning of geometrical size of homogeneously 
distributed volumes responsible for melting. 

As a result of fitting we obtain activation energies which are as follows: 

substance methanol ethanol butanol water 
E^, [eV] 1 0.9 0.8 2.5 

Eb, [eV] 3.5 2.9 7.2 Not used 

According to experimental data (Fig. |2]), water doesn't indicate any significant influence 
of premelting on positronium lifetime especially before its melting point, therefore for it we 
assumed A?, = 0. 
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